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ABSTRACT 

The Damped-Lya systems (DLAs), seen in absorption in the spectrum of back- 
ground quasars, are believed to contain a large fraction of the neutral gas in the 
Universe. Paradoxically, these systems are more difficult to observe at z a b s < 1.7, 
since they are rare and their H I feature then falls in UV spectra, thus requiring the 
use of space-borne facilities. In order to overcome this observational difficulty, Rao & 
Turnshek (2000) pioneered a method based on Mgll-selected DLAs, that is absorbers 
discovered thanks to our knowledge of their Mgll feature in optical spectra. In the 
present work, we use new observations undertaken at the TNG as well as a careful 
literature & archival search to build samples of low redshift absorbers classified ac- 
cording to the technique used for their discovery. We successfully recover JVhi and 
equivalent widths of Fell 2600, Mgll 2796, Mgll 2803 and Mgl 2852 for a sample of 
36 absorbers, 21 of which are Mgll-selected. We find that the Mgll-selected sample 
contains a marginally larger fraction of absorbers with log(A r Hi)>21.0 than seen oth- 
erwise at low redshift. If confirmed, this property will in turn affect estimates of fini 
which is dominated by the highest H I column densities. We investigate on the source of 
the potential discrepancy and find that log(iVHi) does not correlate significantly with 
metal equivalent widths. Similarly, we find no evidence that gravitational lensing, the 
fraction of associated systems or redshift evolution affect the absorber samples in a 
different way. We conclude that the hint of discrepancies in ./Vhi distributions between 
the Mg II-selected DLAs and the others most likely arises from small number statis- 
tics. Therefore, further observations based on both HI and Mgll selection techniques 
are required to better clarify the impact on estimates of J7hi at low redshift. 

Key words: cosmology: observations - large scale structure of Universe - galaxies: 
evolution - galaxies: formation - intergalactic medium - quasars: absorption lines. 



1 INTRODUCTION 

Intervening absorption systems in the line of sight towards 
distant quasars allow for direct observation of the distri- 
bution of gaseous matter from the epoch of initial galaxy 
formation to the present day (Lanzetta et al. 1991; Wolfe 
et al. 1995). In particular, these absorbers provide impor- 
tant ways of measuring the neutral gas and metallicity 
content of the Universe at high-redshifts (Vladilo et al. 
2000; Savaglio 2001; Kulkarni & Fall 2002; Prochaska et al. 
2003; Peroux et al. 2003a). The high-column density end 

* Marie Curie Fellow 

f e-mail: cperoux@eso.org 



of this population is composed of the Damped-Lycv systems 
(DLAs) which have HI column densities, log(A?Hi), >20.3 
cm -2 and the sub-Damped-Lya systems (sub-DLAs) with 
log 19.0<N(HI)<20.3 cm -2 . These two are found to contain 
a large fraction of the neutral gas in the Universe (Peroux 
et al. 2003b). Therefore, these systems can be used to com- 
pute the ratio of HI density to the critical density of the 
Universe, which in turn provides information on the gas 
consumption and star formation over time. 

This cosmological mass density is rather well con- 
strained at high-redshift (Storrie-Lombardi et al. 1996; 
Storrie-Lombardi & Wolfe 2000; Peroux et al. 2003b). In 
contrast, measurements of Ohi at z < 1.7 are paradoxi- 
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Table 1. Journal of the TNG observations for metal equivalent width measurements. 
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a Quasar with Broad Absorption Lines (BAL). 



cally more difficult for several reasons: the observed absorber 
wavelengths are shifted to the ultraviolet requiring space ob- 
servations and the geometry of the Universe combined with 
the paucity of DLA systems requires the observations of 
many quasar lines of sight. Nevertheless, this redshift range 
is of particular importance since it comprises ~70% of the 
look-back time (Om = 0.3; Oa = 0.7 cosmology). 

Lanzetta et al. (1995) were the first to derive fim at low 
redshift using the International Ultraviolet Explorer (IUE) 
satellite. Rao, Turnshek & Briggs (1995) proposed a new 
method based on the observational evidence which indicates 
that DLAs are always associated with a Mgll system, while 
the reverse is not true. From successful HST observations of 
a sample of Mg Il-selected systems (Rao & Turnshek 2000) , 
they estimate DLAs statistics by correcting for the obser- 
vationally known incidence of Mgll systems in a random 
quasar sample. They then compute the mass of neutral gas 
using their "derived" sample of low redshift DLAs. Using the 
same technique, Churchill (2001) has also discovered new 
DLAs at very low redshifts: z ~ 0.05. These observations 
imply a rather flat distribution of the cosmological mass 
density of neutral hydrogen down to low redshift, at odds 
with current ideas and models of cosmic star formation (Pei 
& Fall 1995; Nagamine, Springel & Hernquist 2003). More- 
over, local estimates of H I mass measured by Zwaan et al. 
(1997), and recently confirmed by Zwaan et al. (2003) are 
difficult to reconcile with such high values of f2m in low 
redshift quasar absorbers. 

In the present work, we compare the observed properties 
of Mg Il-selected absorbers with more "classical" DLAs/sub- 
DLAs. The goal is to identify whether both observational 
methods sample the column density parameter space and 
do not impact the calculation of the cosmological mass den- 
sity at low redshift. In section 2, we present the data samples 
which are from both a targeted observational program and 
a dedicated literature and archival search. We describe the 
data reduction process as well as the method used to make 
the measurements of log(A r Hi) and metal equivalent widths 
(hereafter EW). Constructing well-defined samples is chal- 
lenging since the data have been taken at various epochs 



by different authors having distinct scientific goals in mind. 
Nevertheless, we have been extremely conservative in se- 
lecting the systems and, for the purpose of our analysis, we 
have also split the non-Mg II selected samples into two sub- 
samples which allow for comparison of the various discovery 
techniques. Details on this matter are given in sections 3 
and 4 where we also compare the properties of the samples 
of classical and Mg Il-selected quasar absorbers. 



2 THE DATA 

2.1 Mg II Measurements 

In order to build the so-called "classical" sample, we have 
acquired optical spectra covering the Fell 2600, Mgll 2796, 
Mgll 2803 and Mgl 2852 features of known DLAs/sub- 
DLAs with well constrained A^hi column densities. The ob- 
servations were carried out in service mode with the DO- 
LORES spectrograph at the 3.58 m Telescopio Nazionale 
Galileo (TNG) located at Roque de Los Muchachos Obser- 
vatory in the Canary Islands. High signal-to-noise optical 
spectrophotometry was obtained covering the appropriate 
wavelength range so as to study the Mgll and other asso- 
ciated metal lines in z a b s < 1.7 quasar absorbers, the exact 
range depending on the grism used for the observations. A 
journal of the observations is presented in Table 1. 

Thirteen quasars were observed between February and 
April 2003. The integrating times are all of a total of 1800 
seconds regardless of the brightness of the object. Only in 
one case (Q1354+3139/Q1351+318), the weather conditions 
were poor and the resulting spectrum thus suffers low signal- 
to-noise ratio. A combination of high and low-resolution 
modes of the DOLORES spectrograph is used. The cam- 
era is equipped with a Loral thinned and back-illuminated 
CCD with 2048 x 2048, 15 [i pixels. Grisms 1, 5, 6 and 7 
(LR-B, HR-B, HR-V and HR-R respectively) are used 
for the observations. Therefore the dispersion of the spectra 
ranges from 0.8 to 2.8 A pixel -1 . All the observations were 
taken with a slit width of 1.0 - 1.5 arcsec. 
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Figure 1. Fell, Mgll and Mgl measured as part of our survey. Thirteen of these systems were observed at the TNG while two 
additional absorptions were measured in an archival HST spectrum (Q0238+1636/0235+164) and in the commissioning ESO/UVES 
spectrum, degraded to the TNG resolution (Q1106-1821/HE1104-1805). 
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The data reduction is done using the IRAF 1 software 
package. The bias frames are so similar over successive nights 
that a master 'zero' frame for each month is created using 
the IMCOMBINE routine. After trimming the data, they 
are zero corrected using CCDPROC. Similarly a single flat- 
field frame is produced by taking the median of the flats. 
The overall background variation across this frame is re- 
moved using IMSURFIT to produce an image to correct 
for the pixel-to-pixel sensitivity variation of the data. The 
task APALL is used to extract 1-D multi-spectra from the 
2-D frames of the quasars and spectrophotometric standard 
stars, while a cut through the 2-D frames of the Arc served 
as an extraction procedure. The APALL routine estimates 
the sky level by model fitting over specified regions on either 
side of the spectrum. In the case where the object is rather 
faint on the 2-D frame (i.e. for Q0238+1636/0235+164), dif- 
ferent exposures are IMCOMBlNEd before extraction with 
APALL. The spectra are then wavelength calibrated using 
Ar, Ne or He arcs according to the grism used. We pay par- 
ticular attention to this step of the data processing, since 
it is crucial for the later measurements of the metal equiva- 
lent widths. Given that the objects have magnitudes ranging 
from 15 to 19 but identical exposure time, the signal-to-noise 
ratios of the resulting spectra vary from one object to an- 
other. 

In addition to these new data, we thoroughly checked 
published literature and public data sets for any other 
spectra of sufficient quality for us to estimate equiva- 
lent widths of metal features. Fell 2600, Mgll 2796, 
Mgll 2803 and Mgl 2852 lines found in such a way are 
only kept if they are of comparable quality to the one 
obtained with our own TNG survey. Only two spectra 
met those criteria: the UVES commissioning spectrum of 
Q1106-1821/HE1104-1805 and the HST/STIS spectra of 
Q0238+ 1636/0235+ 164 taken as part of program 7294 (PI: 
Cohen). In order to obtain a uniform set of data we de- 
grade the UVES spectrum of Q1106-1821/HE1104-1805 
to the TNG data resolution before undertaking the equiva- 
lent width measurements. 

The Fell 2600, Mgll 2796, Mgll 2803 and Mgl 2852 
absorption features of all the 15 absorbers are shown in Fig- 
ure We compute the equivalent widths of the metal lines 
for the 15 spectra at our disposal. One object in our TNG 
sample (Q1001+5553/Q0957+561) does not cover the Fell 
2600 line while on several occasions Fell 2600 and/or Mgl 
2852 are not detected at the redshift expected for the ab- 
sorber. It can be seen from Figure 1 that the Mg II doublet 
is resolved in all the spectra used for this study. 

2.2 Ahi Measurements 

In addition to the systems presented in the section above, 
we search the literature for suitable log(A'Hi) measurements. 
In order to gather a significant sample, we look for any 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



absorber down to the sub-DLA definition (log(A r Hi)>19.0) 
and with z a b s < 1.7. We only include JVhi measure- 
ments derived from high signal-to-noise HST (STIS or 
FOS) spectra. We therefore exclude measurements based 
on IUE spectra or Nm estimates from 21cm observations 
or other even more indirect methods. Furthermore, we 
thoroughly check ESO and HST archives for any avail- 
able spectra for which the relevant measurements were not 
available in the literature. Eleven JVhi measurements are 
possible thanks to HST/FOS spectra recovered from the 
HST/ESO archives. The UVES commissioning spectrum 
of Q1106-1821/HE1104-1805 provides a further measure- 
ment of JVhi. In particular, we extend the Rao & Turn- 
shek sample by including 7 sub-DLAs which were MgTI- 
selected but that these authors rejected because they had 
log(A r Hi)<20.3. We also include the two Mg Il-selected ab- 
sorbers mentioned by Churchill (2001) and for which the JVhi 
measurements could be made (Q1429+4747/PG1427+480 
and Q0441-4313/PKS0439-433). 

In order to measure JVhi column densities of the 11 ab- 
sorbers available to use, we normalise the spectra in the re- 
gion of interest using spline functions. We then fit Voigt pro- 
files using the FITLYMAN package within MIDAS (Fontana 
& Ballester 1995). We note that for several HST/FOS spec- 
tra, the absorbing feature does not reach the zero flux level 
in its center, even if the line is clearly saturated. We use 
a convolution of the Voigt profile together with the instru- 
mental line spread function to estimate whether this effect 
is real or not. Only in a few cases do we find it to be an 
artifact and we correct the spectrum by subtracting a few 
percent from the continuum level as is standard practice 
(Boisse et al. 1998; Rao & Turnshek 2000). The typical re- 
sulting error bar in log(JVm) measurements is 0.10 and never 
exceeds 0.15, but does not include error in the continuum 
placement. The resulting fits are shown in Figure [5] for the 
eleven systems, while information on the observing details 
is provided in Table 2. 



3 Mg II EQUIVALENT WIDTHS 
DISTRIBUTIONS 

In this section, we study the distribution of the equivalent 
widths of metals in Mg Il-selected absorbers as compared to 
that of more classical DLAs/sub-DLAs. The latter sample is 
further divided into two sub-samples: on one hand absorbers 
selected by their H I feature and on the other hand systems 
selected at 21cm or in one case, because a galaxy is known 
and later turns out to be a DLA. 



3.1 Samples Definition 

The method developed by Rao & Turnshek (2000) makes 
use of the empirical fact that all DLAs studied are found to 
have a corresponding Mg II metal-like absorption (Rao 1994; 
Rao, Turnshek & Briggs 1995). In their study, the authors 
select a sample of 243 2 abs < 1.65, EW Mgll (2796)> 0.3A 
Mgll absorbers taken from various sources in the literature. 
A fraction of these absorbers (87) already has a spectrum 
as part of the HST archives (48 absorbers) while others are 
observed as part of their own HST observing program (39 
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Figure 2. Voigt profile fits of iVfii retrieved from public archives. Ten of these spectra are recovered from the HST/ESO archive while 
the remaining system is measured in the commissioning ESO/UVES spectrum (Q1106— 1821/HE1104— 1805). 
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Figure 2. continued 
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Table 2. Observing details for the spectra retrieved from public archives and used for Whi measurements. 



7 



Coordinate Alternative z^hs ^obs l°g(^Hl) HST/ESO Inst Gratings Selection 

Name Name prog ID (PI) Method 



Q0424H 


-0204 


PKS0421+019 


1.638 


3207 


19.01±0.15 


6577 (Rao) 


FOS 


270 


Mgll 


Q0427- 


-1302 


PKS0424-13 


1.408 


2927 


19.43±0.15 


6577 (Rao) 


FOS 


270 


Mgll 


Q0427- 


-1302 


PKS0424-13 


1.562 


3114 


19.35±0.10 


6577 (Rao) 


FOS 


270 


Mgll 


Q0826- 


-2230 


PKS0823-22 


0.910 


2322 


19.38±0.10 


6577 (Rao) 


FOS 


270 


Mgll 


Q0938^ 


-4129 


Q0935+417 


1.373 


2885 


20.45±0.10 


6237 (Beaver) 


FOS 


270 


HI 


Q1106- 


-1821 a 


HE1104-1805 


1.662 


3236 


20.75±0.10 


60.A-9022 


UVES 




HI 


Q1325^ 


-6515 


4C65.15 


1.610 


3173 


19.76±0.10 


6577 (Rao) 


FOS 


270 


Mgll 


Q133H 


-4101 


PG1329+412 


1.282 


2774 


19.86±0.10 


6577 (Rao) 


FOS 


270 


Mgll 


Q133H 


-4101 


PG1329+412 


1.601 


3162 


19.33±0.15 


6577 (Rao) 


FOS 


270 


Mgll 


Q1429h 


-4747 


PG1427+480 


0.221 


1484 


19.73±0.10 


6781 (Wills) 


FOS 


130 


Mgll 


Q2131- 


-1207 


Q2128-123/PHL1598 


0.430 


1738 


19.55±0.10 


4581 (Bahcall) 


FOS 


190 


HI 



a UVES commissioning spectrum used for both log(A r ni) and Mgll measurements. 



Mgll doublets). Using these data, they measure whether a 
DLA is present or not at the expected absorbing redshift. 

In the present study, the emphasis is on constructing 
well defined DLA/sub-DLA samples: one where the ab- 
sorbers are Mg II-selected, and another for the remaining 
systems, discovered predominantly, but not only, on their 
HI feature. Therefore the latter also includes absorbers se- 
lected at 21 cm or a known galaxy which is further found 
to produce a DLA signature in a quasar spectrum. More- 
over, we extend the work of Rao & Turnshek (2000) (as well 
as our literature search for "classical" systems) down to the 
sub-DLA definition in order to get a more statistically signif- 
icant sample of absorbers. This is done by measuring the H I 
column densities of Rao & Turnshek's (2000) HST program 
for which Mgll absorbers were detected. 

The resulting Mg II-selected sample in the present 
study is thus larger than the original work from Rao & 
Turnshek (2000). It is composed of 21 systems defined as 
follows: 

• 7 DLAs newly discovered in Rao & Turnshek (2000). 

• 7 sub-DLAs in the same observing program and fitted 
by us. 

• 1 of the 4 very low redshift systems of Churchill 
(2001) for which log(A^Hi) is reported in the litera- 
ture (Q0441-4313/PKS0439-433) and 1 other for which 
log(AT H i) is determined by us (Q1429+4747/PG1427+480). 

• 1 sub-DLA (Q0100+0211/Q0058+019PHL938) se- 
lected by Rao & Turnshek (2000) and 1 other DLA 
found by the same group with their new HST program 
(Q1631+1156/Q1629+120). 

• 2 systems (Q0456+0400/Q0454+039 and 
Q1211 + 1030/Q1209+107) observed during the HST 
program 5351 (PI: Bergeron) based on their metal features 
in optical spectra (see Boisse et al. 1998). 

• 1 DLA (Q1624+2345,Q1622+238/3C336) selected 
thanks to its Mg II signature and observed during HST pro- 
gram 5304 (PI: Steidel). 

It should be noted that most of the quasar spectra with 
DLAs in the HST archives are clearly taken because the ob- 
server expects the absorber to be present. As an example, for 
the system along (Q0304-2211/EX-PKS0302-223), the 
HST proposal abstract (ID: 6224; PI: Burbidge) states that 



"This data set will yield high quality absorption line data 
for studies of the damped Lyman alpha absorption system" . 
Therefore, such absorbers are initially not selected on the ba- 
sis of their Mg II feature and we do not include them in our 
Mg II-selected sample even if they have been "re-discovered" 
in the archive search of Rao & Turnshek (2000) . Rather, we 
define the Mg II-selected sample as all systems which have 
been newly discovered based on their metal signature in op- 
tical quasar spectra. 



The remaining systems form the classical sample 
which can be defined as follows: 

• 10 DLAs/sub-DLAs directly discovered thanks to their 
HI signature. These systems also form the so-called HI- 
selected sample (see Section 4). 

• 4 absorbers selected at 21cm and with Nm 
measured in UV spectra (Q0238+1636/Q0235+164; 
Q0813+4813/Q0809+483/3C196; 

Q1232-0224/PKS1229-021; Q1331+3030/Q1328+307/3C286). 

• 1 system for which the quasar-galaxy asso- 
ciation is known prior to the HST observations 
(Q1215+3309/TON1480). 

Altogether, the number of systems studied here is 36. 
A summary of the resulting sample is shown in Table [3] to- 
gether with the properties of both the quasar itself and its 
absorbers. It should be emphasized once more that we only 
include systems for which both Nm and metal equivalent 
widths can be securely recovered. Hence, while during the 
course of our literature search we have known of as many as 
93 DLAs/sub-DLAs with z a b s < 1-7, our criterion reduces 
the current sample to 36 absorbers, 21 of which are Mg II- 
selected. In particular, we choose not to include the two 
Mgll-selected DLAs (log(Ar m )=21.16 and log(Af H i)=21.41) 
recently reported by Turnshek et al. (2004) since these ini- 
tial results might not be representative of the most recent 
HST UV spectroscopic survey undertaken by their group. 



Table 3. Summary of the properties of the 36 2 a bs < 1-7 absorbers (21 of which are Mg Il-selected) making up our complete sample. 
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a Quasar with Broad Absorption Lines (BAL). 
b Gravitationally lensed quasar. 
References: 

(1) This work; (2) Pcttini ct al. 1999; (3) Chen & Lanzctta 2003; (4) Lcdoux ct al. 2002; (5) Rao & Turnshck 2000; (6) Churchill 2001; (7) Boisse et al. 1998; (8) Boisse et al. 1990; 
(9) Lanzetta 1997; (10) Zuo ct al. 1997; (11) Miller ct al. 1999; (12) Turnshck & Rao 2002; 
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3.2 Metal Equivalent Widths 

We compare the Mgll/Fell EW distribution of Mgll- 
selected and classical systems. A simple picture would have 
Mg II equivalent widths correlating with the column densi- 
ties. Indeed, Churchill et al. (2000) observe that at z < 1.6 
the Mgll equivalent widths are dominated by kinematic 
spreads. They propose that these could be directly corre- 
lated to column density, implying in turn a correlation be- 
tween EW and Am- 

In order to address this issue we plot in Figure [3] Mg II 
2796 (top panel) and Fell 2600 (bottom panel) equivalent 
widths for each of the 36 systems. Upper limits and non- 
detections are plotted at zero. 

In their work, Rao & Turnshek (2000) look for DLAs 
only when Mg II 2796 has a rest equivalent width EW>0.3A. 
This is illustrated in the top panel of Figure 01 by the ver- 
tical dotted line. It can be seen that indeed such a crite- 
rion is appropriate for selecting DLAs with log(7VHi)> 20.3 
cm~ 2 (which is precisely the aim of these authors), while 
any higher Mgll EW threshold will miss a fraction of the 
absorbers. This DLA definition is represented by the hori- 
zontal dotted line. In the present work, we extended both 
classical and Mg Il-selected samples to the sub-DLA defini- 
tion in order to get larger samples. We note that even down 
to lower column densities, the rest equivalent width thresh- 
old set by Rao & Turnshek (2000) is still a good tracer of 
quasar absorbers. 

Nevertheless, it is puzzling to notice that the majority 
of absorbers with log(Am)> 21.0 (horizontal dashed line) 
are Mgll-selected, while fewer of these systems have 20.3 < 
log(iVHi) < 21.0. Conversely, in general classical absorbers 
do not have such high column densities and this regardless 
of the way they have been selected (direct HI selection, 21 
cm or known galaxy). 

In Figure 21 we plot the log(A r Hi) column density as a 
function of metal equivalent widths ratios. The top panel 
shows EW Ms u(2796) /EW Ms u(2803) while the bottom panel 
displays EW Mg //(2796) /EW Fe //( 2 600) • We find again that 
both classical and Mgll-selected absorbers span a similar 
parameter range for both these quantities. 
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Figure 3. Top panel: log(Af H i) as a function of the Mgll 2796 A 
rest equivalent widths. The dashed horizontal line corresponds to 
log(A r Hi)=21.0, while the dotted horizontal line is log(AT H i)=20.3. 
The dotted vertical line is the Mgll rest equivalent width lower 
limit that Rao &; Turnshek (2000) use as a criterion for their DLA 
search. Bottom panel: Same for Fell 2600. 

4 Am DISTRIBUTIONS 

4.1 Comparison of the Two Samples 

In this section, we analyse the normalised cumulative Am 
distributions of the Mgll-selected sample and compare it 
with a purely Hl-selected sample. This distribution directly 
impacts measurements of the cosmological mass density, 
fim. Here, the Mgll-selected sample presented in section 
3.1 is compared with the sample of DLAs/sub-DLAs se- 
lected only on their HI. The latter does not include ab- 
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Figure 4. Top panel: log(-/Vni) as a function of the ratio 
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Figure 5. Normalised cumulative distribution of /Vhi for HI- 
selected and Me II-selected absorbers. 



statistics. In any case, the sheer fraction of high column 
density systems will in turn affect the determination of Qui 
computed directly from the integration of the JVhi distribu- 
tion function. 

Further, we find that by comparing systems with 
log(A?Hi)>20.3 only, none of the Hl-selected DLAs have 
log(7V H i)> 21 (out of 4), while 50% (5 out of 10) of the 
Mg II-selected DLAs have log(Af H i)> 21 (see Figure 0. As- 
suming a poisson distribution, the probability of obtaining 
zero Hl-selected DLAs at log(A?Hi)> 21 when expected 2 
(i.e. half of 4) is 13.5%. Therefore the absence of H I-selected 
log(A?Hi)> 21 systems might also be due to statistical fluc- 
tuation. When all classical DLAs are taken into account: 
only 25% (2 out of 8) have Nm > 21 even though these 
have been selected using different techniques. Already Rao 
& Turnshek (2000) note that the comparison of a sample of 
local gas- rich galaxies studied in 21 cm emission (< z >=0) 
with Mg II-selected low redshift DLAs (< z >=0.5) show 
an unexpectedly high rate of occurrence of very large col- 
umn densities. This could be due to redshift evolution or the 
selection methods or a combination of both. 

Therefore, the Mg II-selected sample has a higher frac- 
tion of log(A?Hi)> 21 systems, but only marginally. Given 
the impact of this on Qhi , we try to identify what could be 
the origin of such a discrepancy by testing various scenarii 
in the following sections. 



function of the ratio EW Mg 7/(2796) /EW Fe //(2600) • 



sorbers selected at 21cm or otherwise but solely H I-selected 
systems. Therefore the total sample is now composed of 31 
absorbers, 21 of which are Mg II-selected. Figure^] presents 
the normalised cumulative JVhi distributions for these. The 
Mg II-selected curve is different from the other at the high 
Nm column density end. In particular, there are systems 
with log(7VHi)>21.0 in the Mg II-selected sample while Hl- 
selected systems are never found to have log(TvHi) > 20.8. 
This cannot be due to a bias in the H I-selected sample, since 
that selection method should in principle favour high JVhi 
column densities. A Kolmogorov-Smirnov (KS) test shows 
that the two distributions are consistent at a 0.56 level of 
probability, but clearly we are dealing with small number 



4.2 Gravitational Lensing 

Another effect possibly at play could be gravitational lensing 
of the background quasar by the absorber itself. This could 
amplify the luminosity of the quasars and therefore make 
them preferentially selected as they would be easier targets 
to observe. In order to check whether gravitational lensing 
affects differently the two samples of absorbers, we plot the 
apparent magnitude of quasars for the samples of classi- 
cal and Mg II-selected DLAs/sub-DLAs. Looking at the top 
panel of Figure |S| we find no evidence that the two pop- 
ulations are composed of quasars with different magnitude 
ranges. There might be hints of a trend where higher column 
density absorbers are found in fainter quasars, probably a 
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Figure 6. Top panel: log(TVjji) as a function of the apparent 
background quasar magnitude. Bottom panel: log(TVHi) as a func- 
tion of the ratio of the physical distance to the absorber over the 
physical distance to the quasar. According to gravitational lens- 
ing theory, the quasar magnification is maximised when this ratio 
is 0.5. 



signature of extinction by dust (Boisse et al. 1998), but this 
result is not currently statistically significant. 

In order to test further the strong lensing hypothesis, we 
perform another test. It is known from gravitational lensing 
theory that the magnification effects are maximised when 
the lensing object is midway between the source and the ob- 
server in physical space. In our given cosmology, the physical 
distances can be expressed as: 

X{z)= [ {l+zf [(l + z) 2 (l + 0.3z)-0.7z(2 + z)y 1/2 dz(l) 
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5x10 4 10 5 

zem— zabs [km/s] 
Figure 7. log(A'Hi) as a function of 2 C m — z a bs in km/s 

In the bottom panel of Figure |S| we plot the ratio of 
the distance to the absorber over the distance to the quasar. 
The configuration which would most favour the occurence 
of gravitational lensing would have this ratio equal to 0.5. 
In our study, the two samples of absorbers do not show any 
significant difference. Therefore, gravitational lensing does 
not appear to be a viable explanation for the observed dis- 
crepancy in log(iVHi) distributions of the classical and Mg II- 
selected DLAs/sub-DLAs. This finding is further supported 
by the results of Le Brun et al. (2000) who have used the 
DLA sample of Le Brun et al. (1997) to show that the ampli- 
fication due to the presence of an absorber along the quasar 
line of sight does not exceed 0.3 mag. 

4.3 Associated Systems and Evolutionary Effects 

Another scenario would have the absorbers associated with 
the quasar rather than being truly intervening systems. In 
order to test this hypothesis, we compute the velocity differ- 
ence between the absorbers and the quasars themselves on 
the basis that the former are actually associated systems. 
The z om — z a b s in km/s can be expressed as: 



(2) 



(■^em Zabs ) X C 

(1 + z em ) 

where c is the speed of light in km/s. The result is shown 
in Figure |7| Interestingly, there seems to be a hint of a cor- 
relation between log(iVHi) and (z cm — z a bs) for the Mgll- 
selected population, while this is not true for the classical 
DLAs/sub-DLAs. In any case, we see no differences between 
the two samples which could possibly explain their apparent 
discrepancy in terms of JVhi distributions. 

Finally, we test for potential evolutionary effects by 
analysing log(AHi) as a function of quasar emission redshift 
(top panel of Figure |8J| and absorber redshift (bottom panel 
of Figure Again, we find no obvious trend, and the two 
populations of quasar absorbers appear to sample the same 
parameter space. 
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Figure 8. Top panel: log(A r ni) as a function of quasar emission 
redshift (2 cm ). Bottom panel: log(TVHi) as a function of absorber 
redshift (2 abs ). 



5 DISCUSSION AND CONCLUSION 



wise at low redshift. This property might in turn affect our 
estimates of Qhi which is dominated by the highest H I col- 
umn densities. 

We investigate the origin of the potential discrepancy 
and find that log(A?Hi) does not correlate in a significant 
way with either metal equivalent widths or ratios of combi- 
nation of equivalent widths, quasar magnitude, quasar emis- 
sion redshift, absorber redshift, velocity distance from the 
absorber to the quasar, or position of the absorber with 
respect to the quasar in physical space (a test for the gravi- 
tational lensing hypothesis). We conclude that the marginal 
discrepancy in Nm distributions between the Mg II-selected 
ones and those discovered with other techniques most likely 
arises from small number statistics. 

We note that regardless of the physical reason for se- 
lecting high-column density systems, these absorbers will di- 
rectly affect calculation of Qui at low redshift. In particular, 
comparison of flm calculations dominated by Mg II-selected 
absorbers at z a bs < 1-7 with high-redshift systems selected 
solely on HI might not be appropriate if the two methods 
are not consistent. Furthermore, if there is a discrepancy be- 
tween classical and Mg II-selected absorption systems and if 
it holds at high-redshift, it would imply that Hl-selected 
surveys for DLAs/sub-DLAs at z>2 have been missing a 
large fraction of the neutral gas. 

Since the Mgll method has proven an extremely effi- 
cient technique to find low redshift absorbers, we anticipate 
that forthcoming observations will help extend and refine 
the current results. Therefore, larger samples are required 
in order to better constrain JIhi at z < 1.7 and enable com- 
parison with high redshift results. 
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In this study, we compare the properties of the com- 
plete samples of Mg II-selected absorbers and more "clas- 
sical" DLAs at z a bs < 1-7. Constructing well-defined sam- 
ples has proven a challenging task given the difficulties to 
build a homogeneous H I-selected sample. Nevertheless, we 
find that the metal EW criterion used by Rao & Turnshek 
(2000) to select DLAs traces well the population of these 
systems down to the log(A r Hi)>19.0. Moreover, we find that 
the Mg II-selected sample presented has a marginally larger 
fraction of absorbers with log(7VHi)>21.0 than seen other- 
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